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Inspired 


One day I came accross this formula of Ramanujan (notebooks) 


(3) = Tit” 25" 1 
S(3) = 759 nai n3(e2™ — 1) 


I tried to understand this formula, 

for doing it I made a series of experiments 

with one of my favorite program : lindep or PSLQ, 

that is Integer Relation algorithm. Lindep is part of PariGP 

and now PSLQ is no longer a cryptic FORTRAN animal but part of Maple. 
I made an interface within maple to write a fortran source from 

1 inquiry, compile it on the host computer, run it and come back with the 
answer. I made one for Mathematica and Maple too. 


This lead to more findings 
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You see the pattern here ?, e”” e277” e47™ 


And some exotic ones like 


e” —7 = 19,99909997919 ... in 1987 with my hp 15c and is still a mystery. 


ae 3 


nN 
» —“___ = 10,000000000000000190161767888663 ... 
e mn/7 _ 1 
n=1 
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Pe n 
» —_——. = 119.0000000000000000000000000000000959374585 ... 


ae e2mn/13 = 1 


Getting back to Ramanujan again, an identity with 1. 


a nis 
24) = 7 =i 
n=1 


In fact, there are more like that 


More formulas 


— ee 
24 » e2mn 71 =a 1077 digit prime 
n=1 


And 
= 22807 
240 » pum 4 a 71399 digit prime 
n=1 


Following the lead we get : 


aul oe Fs 
= 16 
eH= 16) aq 7 2 eee 
n=1 n=1 


Can we get other primes like that ? ... all the primes ? 


But why 691 ? 


ue pill a ial 
=. 16 
691= 16 ) m1 2) 
n=1 n=1 


In fact, it comes from this identity with Eisenstein series 
(Jean-Pierre Serre, cours d’arithmétique, p 157.) 
65520. 

jor ZL, eka 


n=1 


Els, — 1+ 


Hum, 2° — 16 = 65520, and 691 is the numerator of 
B,>2,a Bernoulli number. 


Depending of the person, you could prefer the first 
version or the 2". 
Here Eis, is the Eisenstein series (not Euler numbers) 


Anyhow, we have this approximation of 691... 


2711! 
Well, yes there are others 
2°6! 
61 = 


Tt’ 


And 61 is the 3rd Euler number. 


2°12510! 
E510 * tT 


Es19 isa 1062 digits prime, 


These numbers come from the expansion of the Dirichlet 
beta series like 


oo co Ey 
B(s) = ye (2n+1)5 


Where are the Euler numbers coming from ? 


x? 5x* 61x4 
=1+24— + 


cos(x) _ 2! 6! 
re x? 7 5x* 61x4 
cosh(x) 2! 6! 


Mote generally, can we find all trisonometric expressions 
that would lead to approximations of primes with 2* ? 


Could it be possible to get all the primes with these intriguing 
expressions with n* ? 


For example, with another expression (trig) we have 


29 


13 (4 + 22) 


If we look at A006873 (Number of alternating 4-signed 
permutations) 
1, 1, 7,47, 497, 6241, 95767, 1704527, 34741217, ... 
If the sequence contains a prime then if we have the asymptotic 
expansion of a(n) it leads to one more approximation of that prime 


using 7 


7 = 


In this case, itis 


nm" |(442y2) 


The E.G.F. of A006873 is =O *=sG#) 


cos(4x) 


But, how the expression of a(n) is found ? 


1) expand Se ON) te a series 
cos(4x) 

2) Collect coefficients of the e.g.f (with nl) 

3) Compute the ratio of a(n + 1)/a(n) 

4) Compute first differences 


5) Identify the constant : 2.546479089470325372302... 
. 8 
6) The constant is Z 


7) Retro engineer the expression to a(n). 


8) |(442v2)1 is found in ee = |(4+2v2) 


Building one by one each prime from this idea. 


768 64/2 
11 = ae bp 5 
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1944/3 2 ial 
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4-10! 4-8! 
oh Sa ae 
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Not much of a pattern found here 


We have some primes with Euler numbers via the Beta 
Dirichlet series, some Bernoulli numbers... 


Can't we just generate primes with these informations ? 
Not exactly. 


Can't we just generate some primes with <any> formula 
? 


What are the known formulas anyway ? Which is the most 
efficient ? 


Eratosthenes 


Mersenne 


Fermat 


Euler 


Mills and Wright 


Wilson 


Jones, Sato, Wada, 
Wiens 


John H. Conway 


Dress, Landreau 


Benoit Perichon and 
al. 


Tomas Oliveira e 
Silva et al. 


Fridman et al. 


A064648 


Simon Plouffe 


-276 to -194 


1536 


1640 


1772 


1947 and 1951 


Circa 1780 
1976 


1987 
2010 


2010 


2019 


2019 


2019 


2019 


Sieve 


Primes of the form 
2P-1 


Fermat’s little 
theorem 


Second degree 
polynomial 


Double exponential 


Uses p! 


25th degree 
polynomial with 26 
variables 


FRACTRAN 


6th degree 
polynomial 


26 primes in 
arithmetic 
progression 


Sieve optimized, 
fastest known 
prime gererating 
program 


Prime generating 
constant 


Engel expansion of 
0.705230171... 


Efficient Mills- 
Wright-like formula 


Practical 
Practical, exact 
Weak Probable 
primes 
Practical 


Practical 


Theoretical 


Theoretical 


Theoretical 


Practical 


Practical 


Practical 


Practical, limited to 
precision 


Practical, limited to 
precision 


Practical, limited to 
precision 


Computable infinity 


51 

Computable infinity 
40 

Less than 5 known 
exactly 


Very few primes 


Very few primes 


Very few primes 


58 


26 


Computable infinity 


Computable infinity 


Computable infinity 


Computable infinity 


The 6'th degree polynomial took months to find 


Prime numbers are hard to generate. 
The formula of Mills is a good example 


if A = 1.3063778838630806904686144926... 
then |A®" | is always prime. 


a(n) = 2, 11, 1361, 2521008887, 


16022236204009818131831320183, Woops , 
41131011492151048000305295 ee a €ac 
37915953170486139623539759 iteration 


933135949994882770404074832568499, ... 
a(21) is 1.214 billion digits long 


Here is the algorithm of Mills seen in the 
eye of reverse-engineering. 


Begin with p=2 


a) New prime = Next Prime(p‘3) 
b) Go toa). 


You get the sequence, 2, 11, 1361, 2521008887, 
16022236204009818131831320183, ... 


The formula of Wright is even worse 


E. M. Wright formula (1951) 


if do = a = 1.9287800... and g,41 = 29" then 
L9n| = ae | is always prime. 


a(n) = 3, 13, 16381, 

The fourth term is 4932 digits long. 

No one was able to compute the 5" term of 

this sequence. 

In both cases, it is a good idea but not practical at all. 


But if we used Sylvester's sequence rather 


A000058 in the OEIS catalogue is 
2,3, 7,43, 1807, 3263443, 10650056950807, 


has the property that 


1=5+- +- = + —+—— +. 
7 1807 


Called the Pierce expansion or egyptian 
fractions expansion of 0.99999999999... 


The sequence is given by the recurrence 
Sol = Ss =e 
So starting at 2 we get 3, then 7, 43, ... 
But what if we start differently by having 
So = 1.6181418093242092... 
Will produce 2, 3, 7, 43, 1811, 3277913, ... 


all primes. Nice, but it grows too fast. 
The length doubles at each step. 


The number 
So = 1.6181418093242092... 
Was found using simulated annealing + Monte-Carlo 


Simulated annealing is what we call 
« le recuit simulé » 


Simulated annealing + Monte-Carlo 


1) First we choose a starting value and exponent 
(preferably a rational fraction for technical reasons). 


2) Use Monte-Carlo method with the Simulated Annealing, 
in plain english we keep only the values that show primes 
and ignore the rest. Once we have a series of 4-5 primes 
we are ready for the next step. 


3) We use a formula for forward calculation 
and backward. 


One example 


Hypothesis : there exist a infinite sequence of 
primes generated by {c: n"}, c real and { } is the 
nearest integer. 


Yes, if c = 0.2655883729431433908... then the 
sequence (n > 3). 


7 , 67, 829, 12391, 218723, 4455833, 102894377, ... 


But fails after 19 terms at n=22. The sequence is 
finite. 


We go back to Mills model 


What if we use a smaller exponent and test if 
it works ? 


When ap= 43.80468771580293481... then if 
5 


An+1 = QAy*,and use { } to isolate primes. 
This is now sequence : A323176 

113, 367, 1607, 10177, 102217, 1827697, 
67201679, 


Now, if you want an even smaller exponent choosing carefully 
QA, would it work too ? 


di 
Lets try: Q,44 = Qy10 


Then if aj = 
100000000000000000000000000000049.3122107477645 


We have this prime sequence : 
100000000000000000000000000000049 
158489319246111348520210137339236753 
524807460249772597364312157022725894401 
3908408957924020300919472370957356345933709 
70990461585528724931289825118059422005340095813 


Continuing according to this idea... 


101 


Then if ay = 10°°° + 961.49937633785074906096890050... 


I could compute 100 terms of this sequence : a(100) is a 1340 
digits prime (only). 


http://plouffe.fr/Record%20100%20primes%20sequence.txt 


I use a formula for forward calculation and 
backward calculation. 


Forward calculation 
Next smallest prime to { a(n) } 


Backward calculation (to check) 
Previous prime = solve for xin x® — S(n +1). 


This is the simulated annealing 
High speed guessing with a filter 


Guess first value (real value) 


- 0) Apply f(So) 

- 1) Is f(So) prime ? 

- 2) Ifyes, keep the prime in list, if not go to 0 
- with new startng value So. 


- The machine: 
- corei7 at 4.4 Ghz with a 220 TB jbod 
- 283 billion f.p.o. per second 


Finally what could be the sequence with the 
smallest initial value, like 2 ? 


Let’s try :a(n) = {29} 
Where d = 1.3007687041481769105525256... 


(sequence A306317) 


Z, 3, 9, 7, 13, 29, 79, 293, 1619, 1494 75eZGoeer 
11570443, 1540936027, ... 


Can we go backward too ? 
Like from any specific prime number 
using this algorithm ? 


Let’s say from 109° + 267 to 2? 


Yes, if the exponent a is inverted, 

Whena = 0.38562256415290... 
Then we have the sequence : 
10799 + 267, 742123524365563, 542489, 
Gs, /, 2. 


Here a(0) = 2.1322219996628413452 and 
the exponent 
1/ a = 2.5932092490404286167308... 


In 1902, a certain M. Cipolla published a formula 
for the n'th prime number. 


Pn = nUin(n) + Indin(n)) - 1 + o(n) 


On the other hand the formula for the number of 
primes less or equal to n is 


One formula being the functional inverse of each 
other. 


Actually, no. 
Very recently|[4,5], a number of people began to 
realize that these inverses are not as they appear. 


If m(n) = aa ay 
then the inverse is -y W(- 1/y) (or for n) to 
simplify the notation. 


This means that p, = —nW_, (— -). 
Knowing that the value of W(—1/n) has to be 
with W_, and not W (0, — =) 


Now there is a big question about p,, 
m(n) and the precision. 


As we know the PN.T. is a major item. 
But, in term of precision : it is very rough. 


It is true yes, but when n- oo, The same 
with pp. 


We go back to the classic equations. 


n 
1(n) ~ In(n) 


Is the classic equation, we change it for, 
(see Dusart thesis [2010] for details) 


1(n) ~ In(n) —1 


If the compute the inverse 
n 


In(n) — 1 mica 
solve(%,n); gives (y is renamed). 


€ 
Pn ~ —nW_, (= =) 


With the nth prime we have the formula: 


Pn ~ nlog(n) 


Then 
n 
n(n) = ——— 
Wo(n) 
Now we will look at what the error looks 
like for 
Pn ~ —nW_, (—-) 


Then, one by one we eliminate different 
hypothesis about the difference between 
our calculated p, and the real value. 


- The size of what is left is comparable to 
one of the <straight lines> of W;,(n). 

- At first sight the value of the difference 
is a Straight line (correlation is 0.9999), 
It is not. 

- What is significant is the magnitude of 
the difference only. 

- The W,(n) and logarithmic fit are 
indistinguishable. 


For this comparison we need to consider the 
extent of the known tables of p,. 


My own table is up to 20000 billions. 
The known long range table is only up to 1 
(powers of 10 only). 


02+, 


The table for z(n) is up to 27000 billions and 
the long range goes up to 107”. 


The ponctual very large known primes are 
useless for this study for a very simple reason : 
We do not know the rank of these primes. 


n 


for large values of n 


W k(n) 


n 


m(n) is the straight line in blue between the values of ——— 
|W(n)| 


Since the values lies in between, it is natural to use 
differences to realize we still have the same thing but smaller 


nih me on eS 


nN nN nN 
Wats) ~~ teeny Wate) 


This is strange. 


The known behaviour of p,, and (n) is not 
exactly predictable, the evaluation of Riemann 
with li(x) not simple. 


The last computations of p,, and m(n) were 
‘difficult’. (Months of computer time). 


From numerical evidence then 


m =n, (-£)- ps 


The surprise is that : what is left (again) is 
something that resembles exactly what we had in 
the first place (!). The ‘curve’ is still a ‘straight line’ 
but the magnitude is smaller. 


The only plausible explanation is that we have here 
the matryoshka principle: russian puppets. 


Dn © nlog(n) 


piles & = m eS 


is more precise 
Tested at n=10** 


P1024 © 58308642550474983476717666 


The real value being 
58310039994836584070534263 


Now if we continue with this matryoshka principle, 
what is the next term | 


For P4924 by using a bisection method to find the next 
term in the form of 


n 
W,(n) 


The next terms are 114 and 96606. 
W,.(n) can be approximated with formula 4.20 in 
Corless et al. 

W,.(n) = log(n) + 2mki + log(log(n) + 27tki) 


Pio2 = —1074W_, ( : ) a. 


~ 1024) ~ w,(4024) * 


n n 
a + $$$—eE7o + eels 
W114(10**) — Wo¢606(10**) 


The value is 

58310039994824799949493554 compared to 
58310039994836584070534263 

(12 exact digits). 

With 3 terms : 6 exact digits. 


In 1994, B. Salvy published a paper to dig out an algorithm to 
get dozens of terms in the Cipolla formula: 


In(in(n)) — 2 _ In(In(n))? — 6 In(In(n)) + 11 A -) 


Dn~n (ine + In(In(n)) — 1+ In(n) 2 In(n)2 


The formula is quite similar to the asymptotic expansion of 
Wo(n) 


(2 I, L,>(6—9L, +212) L,(—12 + 36h, 22a 
W,(n) es oe 2( 5 ue 2( - a 2( 2 ; 2 De. 
iB 218 abe 127, 


Here L; = In(n) and L, = In(In(n)). 


In principle, with enough terms and withn >> 1 
it should do the thing. 


Not exactly, with 72 terms of Cipolla-Salvy 
formula we get 12 digits exact too. 


There is a limit to it. 


The expansion in Lambert functions is much 
simpler. 


For the moment, the only clue I have about the 4" 
term is that it is proportional to the 
[logz(n) + 1], that is the log of n in base 2. 


For the info, here is the 10" term of the Cipolla-Salvy 
expansion. 


k*In(k)*(1+ (In(In(k))-1)/In(k) + (n(In(k))-2)/In(k) *2+ (-1/2* 
In(In(k))*2+3*In(In(k))-11/2)/In(k)*3+(1/3*In(In(k))*3-7/2* 
In(in(k))*2+14*In(In(k))-131/6)/In(k)*4+(-1/4*In(In(k))*4+ 
23/6*In(In(k))*3-49/2*In(In(k))*2+159/2*In(In(k))-1333/12)/ 
In(k)*5+(1/5*In(in(k)) *5-49/12*In(In(k)) *44+73/2*In(In(k))*3 

é -367/2*In(In(k))*2+3143/6*ln(In(k))-13589/20)/In(k) *6+(-1/6 
*In(In(k))*6+257/60*In(In(k))*5-1193/24*ln(In(k))*4+1027/3* 
In(In(k))*3-17917/12*In(In(k)) *2+47053/12*In(In(k))-193223/ 
40)/In(k)*7+(1/7*In(In(k)) * 7-89/20*In(In(k))*6+959/15*ln(In 
(k))*5-13517/24*In(in(k)) *4+6657/2*In(In(k))*3-39769/3*ln( 
In(k))*2+493568/15*In(In(k))-32832199/840)/In(k)*8+(-1/8*ln 
(In(k))*8+643/140*In(In(k))*7-14227/180*In(In(k))*6+34097/ 
40*In(In(k))*5-76657/12*In(In(k))*4+616679/18*In(In(k))*3-\ 
642111/5*In(In(k))*2+36780743/120*ln(In(k))-893591051/2520) 
/In(k)*9+(1/9*ln(In(k)) *9-1321/280*ln(In(k))*8+119603/1260* 
In(In(k))*7-218809/180*ln(In(k))*6+1328803/120*In(In(k))*5-\ 
2696687/36*In(In(k)) *4+33904723/90*ln(In(k)) *3-40633409/30* 
In(In(k))*2+7921124011/2520*ln(In(k))-2995314311/840)/In(k) 
“10) 


But, let’s go back to m(n), we had 


We apply the same scheme let’s say for 10” for z(n) and 
Pn - 


D107 —e 1 1 1 i 
107 © Gay) ~ Wo (107) ty, (107) ty. (107) ty, (107 
0) 22 763 5323546 


meu) 1 1 A 1 ri il 
107) -W,(107) W,.(107) W.4,007) wees 


It is similar, how similar is it ? 


Graph of the coefficients of m(n) and p,, in the 
Lambert expansion, every billion from 10? to 19674 : 10”, log scale. 


3 terms of p(n) and pi(n) , Lambert function expansion, log scale 
from 10%9 to 19674*10%9 
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—_— p(n) is term —— pin) 2nd term —— pin) 3rd term —— pi(n) Ist term —— pi(n) 2ndterm —— pi(n) 3rdterm 


This is where we can apply the duck principle. 


If an animal has a beak like a duck , feathers like a duck, the 
color of a duck, quacks like a duck and has 2 feet like a duck 
then it’s a duck. 


With this Lambert function expansion the 2 quantities z(n) 
and py, are the same (except for the first term of p,,). 


Recently, a certain André LeClair and Guillerme Frang¢a (2014) 
had a formula for the nth zero of Riemann’s Zeta function. 


It follows the same idea. If N(n) is the number of non-trivial 
zeros (considering only the imaginary part) then 


\ 11 


N(n) & = log Oe — 


— >See 
2 ac 
By inverting (functionally) the formula we obtain a formula for 
the nth zero. 
(8n —11)z 


The formula is spectacular in precision. 


3 TU 
o(1)¥ —~ ——,~ 
* w (ge) 
Se 
Is 14.5213469... when the real value is 14.13472514 


So precise that they could evaluate precise values of 
o(n) with n = 101° by using an additional Newton- 
like interpolation. 


We have here a quantum leap compared to previous 
models. 


Again, if we go back to the classic known equations. 


And N(n) (Riemann) is 
n n 
N(n) ~ = log() 


27e 


And now by solving for n in each case we get 


N(n) ed a(n) | 7 —20 
2mn w(Z) 2mn -( 
e 


Pn e 1 1 
Se —W_ so 
n : ( = Wor) Weoto mj+i™) 


If we collect the 4 formulas we found, dividing by 
either n or 27n we get 


a 1 1 > i” —e 1 
m(n) = "(eam Pn =n(- a()-a5-) 


= 2 W. mae N(n) = 2mn — 
On — 47N (- -( 5 ) a (®) W,(n) 4 


From there, 2 possibles directions 


If the Euler principle applies then we should 
have a sum and a product on each side. 


Or, the expression with primes needs to be completed 
with an expression using the zeros of the Zeta function 
(1/2 + it), then it has to match with the equivalent 
expression with the primes on the other 2 equations. 


I leave this question as an exercice to be completed... 


Here is a model we can try for z(n) and p, 


mn) 1 1 
n ~ Wo Wo(n) — = Mil (n) + ekf@™ 
—e 1 
= —_w_,(—)—-——_-- ) —=_—ame 
Pn ( n TAGS Lu W,.(n) + ekg 


Here f(n) and g(n) are found using a logarithmic fit. 


For f(n), at n = 4635000018752, the formula gives 
146388867645773 exactly. 


For g(n) , the formula is quite similar. 


Comparing the 2: 


f(n) = 5.1407131338852538860618655508885 + 
.048089483129908800105508959416636 In(n) 
g(n) 

= 5.1425259035418911897661770856362 

+ .047839197978255729085511308229899 In(n) 


Is this just a coincidence ? 


Let’s try another model ? 


mn) 1 1 
n ~ Wo Wo(n) = Wy(n) + kim 


~ —W (—) : bi 1 
ee An W,(n) £4 W;,(n) + k9™ 


g(n) = 4.8493349460 + 0.0557287326 In(n) 
f(n) = 4.8604042576 + 0.068197411 In(n) 


I think that there is a possibility that one could find 
eventually an exact formula for m(n), p, , N(n) and 


a(n). 


And now, something completely different. 


This is Viéte formula (1593). 


2 2 2 2 
Ja Jzaws l2+Voeva J2+\2+v2+v2 
2 2 2 2 

As far as mathematic is concerned, nobody knows what is 


2 — 
the binary expansion of . This is just a bunch of zeros and ones 


at random. Perhaps we will never find the patterns in it. 


But here is another approach to the problem. 


If we are looking at individual bits of this number 
we do not see anything. 


But let's consider this instead. 

We align all the partial products in Viéte expansion 
and look at the bits as a whole, everything is computed 
in binary. 


The first line is 2 


The 2" line is 2 ane? 


2 escy 2+ 2+v2 


2 


The 3" line is — 


61 


With 99715 bits x 25000 terms gives this : 


Mit. 


An image of 2.5 billion pixels. 
The long vertical lines at the left are the slow 


2 
convergence of the product to - 
We don't see much of a pattern here. 


As a comparison, this is an image of the first 100 millions 
decimal digits of Pi, colorized with 10 colors, blue, green, 
red, yellow... the image is 10000 x 10000. 


If we zoom in: this is pretty much random data. 
One experiment was done with the first 1000 billion digits 
of mz : 10000 images of 10000 x 10000 (see web page). 


I have no explanation for it. 

The successive square roots are producing the effect. 
The thing we can Say is that : from one term to the other 
: bits are not random and the pattern is quite persistant. 


Are there any other algebraic curiosities like that ? 
When experimenting with square roots and square 


roots of square roots, I wanted to know if a persistant 
pattern occurs in the Mandelbrot set. 


As we know, there are similarities. | wanted to know 
what if we approach one limit point ? If the limit 
point is algebraic, are there any patterns in these 
algebraic numbers ? (in binary perhaps ?). 


Well, I have found one formula: 


V167" +1 


16"%4+1 


If f(n) =1+ 


then the binary expansion of f 


has a very, very long and persistant pattern. 


Whenn = 4096 then at position 1342238724 there 
are 4118 successive bits of this number that are all 
‘0’. 


(estimated) at n = 1000000 the persistant pattern 
goes up to the 80000015000004'th position. 
At n= 1000000000 the position is at 8.0 x 107°. 


V 2-48192 4.9,/46819244 


216386 


The first 270 million bits of 1 + 


In false color: 
The bits were colorized to 
enhance the contrast. 


This is the binary expansion of f(100) (same as above) 


f(100)=1.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000010 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000001111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111011000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000010100111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111001010011000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000100101111110111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 
1111111100011010010111100000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
0000000000000000000000000000000000000000000000000000000000000000000000000000000010110101010111110011111111111 
111111111111111111111111111111111111111111111111111111111111111... 
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The first 446 million bits of the number 
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At another point (depending on the choice of width) 


Zoom at the bit level (from an image of 332 m pixels.) 
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The data 


® In all, 74 TB of mathematical data, mostly numbers. 

= 41 TB of primes, from 2 to 80594098476893 (2602 billion entries) 
= 5.6 TB of ¢ zeros, 103 billion zeros 

® OEIS tables (and extended tables) 


= The Inverter, 41 digits (small version with 11.3 billion entries), 64 
digits, 17.2 billion entries 


® Inverter 41 : 1.008 TB, http://plouffe.fr/ip/ 
™@ Inverter 64: 2.15 TB. 
® High resolution images : 1206000, 1.773 TB. 


Errors found... 


® Prime|8200000000]| and Prime|[93000000000] just hangs 
in Mathematica. 


-> The tables of zeros at http://www.lmfdb.org/zeros/zeta/_do contains errors 


when the decimal expansion finishes by 00. 11 errors where found. 


The index and the corresponding zeros do not match. 


In both cases, they are working on the problem. 


Thank you for your attention 


Merci de votre attention 
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